Nonconservative kinetic exchange model of opinion dynamics with randomness and 

bounded confidence 
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The concept of a bounded confidence level is incorporated in a nonconservative kinetic exchange 
model of opinion dynamics model where opinions have continuous values £ [—1,1]. The characteris- 
tics of the unrestricted model, which has one parameter A representing conviction, undergo drastic 
changes with the introduction of bounded confidence parametrised by 8. Three distinct regions are 
identified in the phase diagram in the 8 — A plane and the evidences of a first order phase transition 
for 8 > 0.3 are presented. A neutral state with all opinions equal to zero occurs for A < A C1 ~ 2/3, 
independent of 8, while for A C1 < A < X C2 (8), an ordered region is seen to exist where opinions of 
only one sign prevail. At A C2 (8) , a transition to a disordered state is observed, where individual 
opinions of both signs coexist and move closer to the extreme values (±1) as A is increased. For 
confidence level 8 < 0.3, the ordered phase exists for a narrow range of A only. The line 8 = is 
apparently a line of discontinuity and this limit is discussed in some detail. 
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I. INTRODUCTION 

Many social phenomena can be modelled and analysed 
using the methods of statistical physics. Opinion dy- 
namics is one such phenomenon where the microscopic 
interactions taking place between individuals can lead to 
cooperative behaviour of the society, like the emergence 
of a consensus or agreement. Recently, several opinion 
dynamics models with continuous opinion have been pro- 
posed in which the opinions are updated after a pair of in- 
dividuals interact and exchange opinion. Here the inter- 
action scheme is similar to a kinetic theory of gas, where 
one assumes that gas molecules can collide and conse- 
quently exchange energy and momentum. Two basic 
opinion dynamics models having such kinetic exchange 
scheme are due to Deffuant et al [l] and Lallouache et 
al |2j, which involve completely different concepts. In 
this paper, we attempt to combine the concepts used in 
these two basic models in a single model. In section II, 
we briefly review these two models. In section III, the 
proposed model is described. The results are presented 
and discussed in section IV. Summary and conclusions 
are given in the last section. 



II. A BRIEF REVIEW OF TWO KINETIC 
EXCHANGE MODELS 

Deffuant et al [l[ introduced a simple model (DNAW 
model hereafter) in which opinion exchanges between two 
agents take place only when the difference in the original 
opinions is less than or equal to a preassigned quantity 
5. If Oi (t) is the opinion of the ith agent interacting with 
the jth agent at time t (with \o% — Oj\ < S), then in this 



model the opinions evolve according to: 
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Here 7 is a constant (0 < 7 < 0.5) called the convergence 
parameter and Oi lies in the interval [0,1]. The dynamics 
is such that the opinions tend to come closer after in- 
teraction. Hence as the dynamics proceeds, convergence 
to a finite number of opinions happens; opinions cluster 
around a few values and individuals belonging to different 
clusters no longer interact. The initial distribution of the 
individual opinions is uniform and therefore symmetric. 
This symmetry is broken as the time evolved distribu- 
tion has a multipeaked delta function form. When there 
is only one peak in the final distribution, it is said to 
be a case of consensus, two peaks imply polarisation and 
the existence of a larger (finite) number of peaks signifies 
fragmentation in the society. 

The model is conservative as total opinion is con- 
served in each interaction. Obviously, in this conservative 
model, consensus would imply that opinions converge to 
the value 1/2. Several models have been formulated in- 
corporating the idea of bounded confidence later 
and a general form of kinetic exchange type model pro- 
posed in @. 

More recently, a model in which kinetic exchanges take 
place with randomness, and where there is no conserva- 
tion, has been introduced by Lallouache et al @ (LCCC 
model hereafter). Any two agents can interact in this 
unrestricted model. The opinion evolution here follows 
the rule: 
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where e, e 1 are drawn randomly from uniform distribu- 
tions in [0, 1]. In this model A is a parameter which is 
interpreted as 'conviction'. The opinions are bounded, 
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i.e., — 1 < Oj(t) < 1; in case Oi exceeds 1 or becomes less 
than —1, it is set to 1 and —1 in the respective cases. 

It is possible to rescale the opinions in the DNAW 
model so that they lie in the interval [—1, 1]. Continu- 
ous opinions are relevant in cases like supporting a issue, 
rating a movie etc. Thus setting the interval as [—1, 1] 
appears to be more meaningful since a positive (negative) 
value of the opinion will mean liking (disliking) the mo- 
tion. The magnitude of the opinion would then simply 
correspond to the amount of liking or disliking. For the 
rest of the paper, we thus consider opinions € [—1 : 1]. 

As there is no conservation in the LCCC model, the 
average opinion given by m = \ Y^i °i\/N for a popula- 
tion of N agents, evolves in time and to can play the 
role of an order parameter. It is, in fact, analogous to 
magnetisation in magnetic systems. One can say that 
there is order/disorder when m(t — > oo) converges to 
a nonzero/zero value. The model shows a unique be- 
haviour: below a critical value of A ~ 2/3, all opinions 
identically turn out to be zero while above it, there is a 
nonzero value of the average opinion. Thus for A > 2/3, 
an ordered phase exists. Interestingly, the opinions in the 
ordered phase have either all positive or all negative val- 
ues. Generalisation and variations of the LCCC model 
have been considered in some subsequent works 04§|- 

Symmetry breaking has different connotations in the 
conserved DNAW model and the nonconserved LCCC 
model. In the former, if opinions are initially in the 
interval [—1,1], a consensus implies convergence of all 
opinions to zero value and this is regarded as symmetry 
breaking as mentioned earlier. In LCCC, the identical 
state of all zeros is also obtained below A ~ 2/3 even 
without putting any restriction on the interactions. How- 
ever, this state has been interpreted as a symmetric state 
0. This is following the idea that as to = here, it is 
like a paramagnetic state (which is a symmetric state in 
magnetic systems). But obviously, this is a very special 
paramagnetic state which also has zero fluctuation. 



III. THE MODEL 



ety is fragmented because of lack of communication. A 
large value of conviction is not expected to help reach 
consensus when the confidence level is low when we com- 
bine the two concepts, as the groups stick to their original 
opinions in this case. On the other hand, even when con- 
fidence level is large, very large values of conviction can 
make the society polarised as both positive and negative 
opinions can prevail simultaneously. 

In the model proposed in the present paper, we follow 
eq © for the evolution of opinions but put the restriction 
that agents interact only when |o, — Oj\ < 26. 6 is once 
again the parameter representing the confidence level and 
can vary from zero to 1. 

We therefore have two parameters in the model, A and 
S. 5 = 1 is identical to the original LCCC model. 5 = 
is an interesting limit. Here, agents interact only when 
their opinions are exactly same. We will discuss this limit 
in greater detail later. 

The order parameter in this model is the one defined 
for the LCCC model, i.e., to = \ To avoid 

confusion, we adopt the following terminology for the 
present nonconserved model: when to = and also the 
individual opinion distribution is a delta function peaked 
at zero, we call it a neutral state. If it is not a delta 
function at zero but m = 0, the state will be termed a 
disordered state. When to ^ 0, it is an ordered state; 
further if all individuals have identical opinion (which 
is perhaps only ideological as opinions are continuously 
distributed), it is a consensus state. Hence consensus is 
not merely an agreement in this terminology. Obviously, 
in the conserved system, the nomenclature of order and 
disorder is irrelevant. We intend to locate the regions of 
neutrality, order and disorder in the present model in the 
5 — A plane. 



IV. RESULTS AND NUMERICAL ANALYSIS 

We take a population of N agents and let it evolve 
according to the dynamical rules defined above (i.e., eq 
[5] with a bounded confidence level 6). 



Having conservation in the opinions is rather unre- 
alistic but the concept of having a bounded confidence 
level is relevant in many cases. On the other hand, it is 
true that people may have a conviction level in their and 
in others' opinions. We thus combine the LCCC model 
and the DNAW model by putting the restriction of the 
bounded confidence in the former. It has to be realised 
that the conviction and bounded confidence are indepen- 
dent concepts and the parameters associated with these 
may be regarded as independent. Putting a bounded 
confidence may be applicable to cases in which there are 
groups of people (e.g., members of a political party) who 
share similar opinion and interact within themselves pri- 
marily. On the other hand, conviction is associated with 
individual behaviour. It is known that in models with 
bounded confidence, for small confidence level, the soci- 



A. Order parameter 

The behaviour of the order parameter after the system 
reaches equilibrium is presented in Fig. [T] As a func- 
tion of A, we find that the order parameter first assumes 
non-zero values at a threshold value of A = A Cl which 
is independent of 5; A Cl ~ 2/3 as in the original LCCC 
model. The order parameter increases with A beyond A Cl 
up to a certain value of A and decreases to zero as A is in- 
creased further. The decrease becomes steeper with d and 
more so when the system size is increased. The results 
indicate that there are three distinct regions: one ordered 
region for intermediate values of A and two regions at low 
and high values of A where the order parameter vanishes. 
These two regions may be either disordered or neutral. 



FIG. 1: (Color online) Variation of the order parameter m 
with A for 8 = 0.6 for N = 1048 and 4096. Inset shows m and 
s for 8 = 0.2 and 0.4 against A for N = 1048. 

We also study the behaviour of another quantity s de- 
fined as 

s = (\f+-f-\), (3) 

where /+ denotes the fraction of population with opinion 
greater than zero and /_ = 1 — /+ in a particular con- 
figuration. (• • •} denotes average over all configurations. 
It can be easily seen that s is equal to unity both in the 
neutral state and ordered state of LCCC (in the numer- 
ical study, we take an opinion value to be equal to zero 
when its absolute value is less than 10 -8 [Hi). Devia- 
tion of s from unity will indicate that opinions with both 
signs are present in general. We notice that s remains 
close to unity as A is increased from zero before showing 
a sharp fall close to a value of A where the order param- 
eter also starts to fall (Fig. Q] inset). Evidently, as the 
system enters the disordered state, individual opinions 
are > and < with equal probability, s has a mono- 
tonic behaviour and is useful to locate the transition to 
the disordered state. Moreover, comparison of s and m 
shows that these two measures become closer and tend 
to merge as A is increased further. This indicates that 
as one moves deeper into the disordered region, opinions 
become more and more close to the extreme values ±1, 
leading to a polarisation tendency in the opinions. 

B. Individual opinion distribution 

To understand the nature of the phases, the distribu- 
tion of individual opinions may be studied. Such studies 
are known to lead to a correct speculation about phase 
transitions [TTj . This study shows (Fig. |2J that the 
probability for zero opinion is nearly equal to 1 below 
A C1 ~ 2/3 as in the LCCC model for all 6. Hence a neu- 
tral state exists here as well and the confidence level is 
absolutely irrelevant as A Cl is independent of S. 

It may be noted that a plateau like region exists in 
the distribution where the individual opinions have more 
or less the same probability. The extent of the plateau 
depends on A and shrinks with increasing A. One can 



FIG. 2: (Color online) The distribution of individual opin- 
ions shown for N = 2048 different values of A for 8 = 0.4 
averaged over different configurations. Inset shows the same 
for a single realisation of the system. For A < A C1 ~ 2/3, 
the average distribution is a delta function at zero. For 
Aci < A < Ae 2 — 0.7183 (for 8 = 0.4), for a single con- 
figuration, opinions are all of one particular sign, while for 
averages over all configuration, the distribution is symmetric. 
For A > A C2 , the distribution is symmetric even for the single 
configuration. 

investigate the reason for the existence of such a plateau 
region by some further analysis. Let us take specifically 
the example of A = 0.74 and 5 — 0.4 where the plateau 
extends from about \o\ = 0.75 to nearly \o\ = 1. It should 
be noted that even in equilibrium the individual opin- 
ions undergo dynamics. We tag individuals with abso- 
lute value of opinion greater than 0.75 at a time when 
the system has reached equilibrium and note their evolu- 
tion. First, the probability Pt that these remain greater 
than 0.75 in subsequent times is calculated. For a ergodic 
case, i.e., if the opinions can span the entire opinion space 
with equal probability, this probability should be around 
0.25. We find, however, that P t is nearly 0.9 for A = 0.74 
and 6 = 0.4. This means that the opinions remain within 
this interval with a high probability and one need not re- 
ally bother about opinions with lesser values to analyse 
the situation in the plateau region. P t may be close to 
unity if opinions reach a fixed point (or weakly fluctuate 
about one) and does not imply that a plateau like region 
will emerge. Hence we probe further by splitting the in- 
terval into further subdivisions. We estimate separately 
Pi, P 2 and P3 where, 

(i) Pi is the probability that opinions at later times are 
greater than 0.99 

(ii) P2 is the probability that opinions at later times lie 
between 0.89 and 0.99 

(iii) P3 is the probability that opinions at later times lie 
between 0.79 and 0.89. 

One can define other probabilities but these three are 
sufficient for our purpose in the present case. We also 
take three mutually exclusive subsets of the tagged agents 
with initial values of \o\ between 0.75 to 0.80 (subset A), 
between 0.85 to 0.9 (subset B) and between 0.95 to 1.0 
(subset C). Here initial value implies the value at the 
time of tagging after the system has equilibriated. We 
separate out Pi as the probability of individual opinions 
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FIG. 3: (Color online) The probabilities P\,P<2 and P3 are 
shown against A for subsets A, B and C. The data for for the 
different subsets coincide for each of the three probabilities. 
Inset shows the total probability that starting with a opin- 
ion with absolute value greater than 0.75, it remains so at 
subsequent times. Data is for 8 = 0.4 and N = 2048. 



having values very close to unity is much greater than 
the plateau height (this is simply because of the imposed 
boundary condition on the opinions) as shown in Fig. [2] 
The idea is that if opinions do reach (nearly) fixed point 
values, Pi and P3 should be grossly different for all the 
subsets A, B and C (e.g., for A, Pi will be less than P3 
etc.) On the other hand, if P2 and P3 turn out to be 
close to each other for any subset, one can conclude that 
the subspace 0.79 < |o| < 1.0 is uniformly accessible to 
the agents within that subset. Moreover if these proba- 
bilities are independent of subset A, B and C, one can 
argue that a plateau region will exist. We find all the 
three probabilities Pi, Pi and P3 to be independent of 
the subset considered (Fig. [3]). P\ turns out to be much 
higher than the other probabilities as expected and in- 
deed, Pi and P3 are very close to each other for A = 0.74. 
Hence the plateau exists and the reason is the opinions 
here evolve in such a way that only part of the opinion 
space is spanned but spanned with equal probability. We 
have in fact obtained P\,Pi,P-$ and Pt for other values 
of A for 5 = 0.4 and shown them in Fig. [3] 

As A is increased beyond A Cl , in a single configura- 
tion, only all positive or all negative values are obtained 
as N — > 00, while the average over all configurations is 
symmetric about zero as expected. However, as A is in- 
creased further, the opinions, even in a single configura- 
tion, assume both negative and positive values symmet- 
rically (Fig. [2] inset). These results are consistent with 
the results for the quantity s presented earlier. Hence 
we infer that an order-disorder transition is taking place 
at a value A C2 > A Cl which is later confirmed from more 
detailed analysis. For 8 = 1, the LCCC model, A C2 is 
equal to 1 as expected. 



C. Phase transitions and critical properties 

The ordered, disordered and neutral regions may be 
identified in a phase diagram in the 8 — A plane. To 
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FIG. 4: (Color online) The Binder cumulant is shown for 
different values of 8 for N = 512, 1024, 2048 and 4096. For 
smaller values of A, it shows clearly that a negative dip exists 
which approaches the intersection point as N increases. Other 
analysis suggests that the negative dip will be observed for 
all values of A > 0.3 if sufficiently larger system sizes are 
considered (see text). Colour code is same for all the figures. 



obtain the phase boundaries in this plane, we estimate 
the phase transition points by traditional methods, i.e., 
attempt finite size scaling for the relevant physical vari- 
ables, if possible. Among these variables is the fourth 
order Binder cumulant (BC) defined as 



U=l- 



3 (to 



2\2 ' 



(4) 



Here we discuss the case for 6 < 0.3 and 6 > 0.3 sepa- 
rately for reasons which will be clear later. 

6 > 0.3: Plotting U against A, we find that there is 
indeed a crossing point but interestingly, the BC shows a 
negative dip (Fig. [4j) for <5 < 0.6 for the system sizes con- 
sidered. In fact it becomes more negative as the system 
size N is increased and the location of the negative dip 
approaches the crossing point as well. These are typical 
indications of a first order phase transition (T^ |. To con- 
firm whether a first order transition is really taking place, 
we also plot the distribution of the order parameter very 
close to A C2 . One expects the distribution to have peaks 
at nonzero values of to below the critical point (usually 
the distribution is a double gaussian). For a continu- 
ous phase transition, as the critical point is approached 
from below, the peaks occur at smaller and smaller val- 
ues of to, finally merging at to = continuously at the 
critical point. For a first order phase transition, on the 
other hand, the peaks at nonzero values of the order pa- 
rameter remain at more or less the same positions up to 
the transition point [l2Tfl6l |. Here we find exactly this 
behaviour (Fig. [5J; note that for finite systems, weak 
peaks will still show at nonzero values of to just above 
the transition point (instead of a perfect gaussian with 
mean zero). 

We attempt to obtain scaling forms for the BC (U), 
order parameter (to) and a quantity analogous to sus- 
ceptibility per spin (in magnetic systems) given by x = 
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TABLE I: The transition point and the values of the ex- 
ponents for A > 0.3. The typical errors in the data are 
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FIG. 5: (Color online) Distribution of the order parameter 
close to the transition point A C2 — 0.7183 for 8 = 0.4 (N — 
2048). 




FIG. 6: (Color online) The variation of m and x are shown for 
different system sizes for 8 = 0.4. The insets show the scaling 
plots with m and x/ N plotted against x = (A - 0.7183) iV ' 53 . 



jj[{M 2 ) - (A/) 2 ] where M is the total opinion, M 
| Oj | . The expected behaviour are given by 
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(5) 



For first order phase transitions in finite systems, one 
expects that instead of a delta function behaviour at the 
transition point, there will be a peak in the susceptibility 
which will diverge with the system size. The order pa- 
rameter exponent a is expected to be close to zero. We 
find that the above scaling forms are indeed appropriate 



in the present case, the data collapse to a single curve 
for specific values of a, b and /i (shown for m and x m 
Fig. IS]). The value of a is indeed very close to zero and 
b ~ 1 for all values of 5 > 0.3. However, the value of /x 
appears to have a systematic variation with S. Since S 
effectively puts a restriction on the number of compatible 
neighbours, it is not surprising that /i, which is associ- 
ated with N, shows a dependence on S. The values of 
the transition point A C2 and the exponents are presented 
in Table |U 



All the above discussions are however, valid for S > 0.3 
only. The first order phase transition is most strongly 
observed close to 8 — 0.4. As for the negative dip, it is 
not observed for S > 0.6 with N < 4096, but the values of 
the exponents indicate that the transition is first order- 
like. The negative dip for 5 > 0.6 is thus expected to be 
observed for even higher values of N [l2[ ■ 

S < 0.3: When 5 is decreased below 0.3, the results 
do not give any clear indication about the nature of the 
phase transition and shows some anomalous behaviour. 
A rather uncharacteristic behaviour of the order param- 
eter and the Binder cumulant is observed. The order pa- 
rameter m shows a nonmonotonic behaviour when plot- 
ted as a function of 6 with fixed A or vice versa (Fig. [7]) . 
A hump appears in the m versus A plots for S < 0.3 and 
A > 0.7 showing the existence of a local maximum value. 
A closer examination reveals that this hump shrinks in 
size, albeit very slowly, when N is increased for A > 0.71 
(see inset (b) of Fig. [7]) and is expected to disappear 
in the thermodynamic limit. The same is true for the 
hump appearing in the plot of m versus A for S = 0.25. 
For 5 < 0.3, there are large fluctuations and irregulari- 
ties in the BC as well which does not allow one to do a 
finite size scaling analysis and get the exponents. The 
irregular behaviour of the BC and the order parameter 
may be because of the fact that the interactions become 
less likely to occur as the confidence level is decreased. 
For this reason we restrict our study to S > 0.1. 

We have estimated, somewhat approximately, A C2 for 
5 < 0.3 from the crossing point of the order parameter 
curves for different sizes. For example, from the data 
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of m against A for 6 = 0.25 shown in Fig. we get 
A C2 w 0.7. The complete phase diagram in the S — X 
plane is shown in Fig. [8j 



D. Discussions 

Let us try to understand the results by analysing the 
role of the confidence level 5. Let the ith agent with 
opinion Oi interact with another agent with opinion Oi + x 
where \x\ < 28. Then, 
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Oi(t + 1) = A(l + e)oi(t) + Xex. 



(6) 



Consider the case when 8 is small. On an average, 
when A is smaller than 2/3, the first term will make Oi 
smaller in magnitude while the contribution of the second 
term can be neglected. Then we find that Oi(t) — » as 
t — » oo implying the convergence to a neutral state. Since 
for 8 = 1, it is already known that there is a transition to 
the neutral state at A ~ 2/3, we conclude that for any 8 
this is the case as is confirmed by the numerical results. 
Actually, when individual opinions decrease towards zero 
because of the effect of the first term, the difference in 
opinions automatically becomes less (x — > 0) so that large 
and small 8 values have the same effect; the second term 
does not contribute eventually. Thus A C1 is independent 
of 8. 

What happens at higher values of A? First consider 
small values of 8 again which means x is small too. Hence 
the second term still contributes less compared to the 
first and Oi will retain its original sign in most cases if 
A is sufficiently large. So it is expected that there will 
be a region where opinions of both signs are present and 
m = 0, as originally opinions are uniformly distributed 
with positive and negative signs. 

However, if 8 is large, there is no guarantee that the 
second term is small and opinions will retain their original 
signs, unless A is also very high. This explains why we 
observe the transition to the disordered state at a higher 
value of A as 8 is increased. Also, it is not surprising that 
there will be a ordered region between A Cl and A C2 (as 
already it is known to be present for <S = 1.0) where the 
LCCC property of all opinions assuming the same sign is 
still valid. 

Although we have restricted to 8 > 0.1 in the numeri- 
cal simulations, the case when 8 is exactly equal to zero 
can be discussed theoretically to some extent. If an inter- 
action takes place, the opinion for the ith agent follows 
the evolution equation 



0i (t + l) = X(l + e)oi(t). 



(7) 



This equation is nothing but the dynamical equation ob- 
tained for the LCCC model in the limit of a single pa- 
rameter map 0, [13, [TJ[ , where the transition to an or- 
dered state occurs at a value of A = e/4. However, in 
the present model the above is valid only when there is a 
second agent with opinion equal to Oi as well. Since this 
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FIG. 7: (Color online) Variation of the order parameter m 
with A for 8 — 0.25s. Inset (a) shows the variation of m with 
8 for two fixed values of A. Inset (b) shows a magnified portion 
of inset (a) to show the finite size effects. All data are shown 
for four different values of N following the same colour code. 
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FIG. 8: (Color online) The phase diagram in the 8 — A plane 
shows the existence of the neutral region (for A < A c i ~ 2/3), 
the ordered region and the disordered region. The ordered 
and disordered regions are separated by a first order boundary 
(continuous line in red) for 8 > 0.3 obtained using a finite size 
scaling analysis. For 8 < 0.3, the phase boundary (broken 
line in blue) has been obtained approximately only from the 
behaviour of the order parameter (see text). 



will be extremely rare, effectively most of the opinions 
remain frozen and the single parameter map is not rep- 
resentative of all the agents' opinion evolution. In fact, 
8 = may be regarded as a line of discontinuity in the 
phase diagram of the model in the <5 — A plane as it will 
neither have the neutral state nor the ordered state any- 
where. The disordered state is also different in nature 
for 8 = 0; here the individual opinion distribution will be 
flat while for 8 0, however small, it is not so. 



V. SUMMARY AND CONCLUSIONS 

In summary, we have studied a model of continuous 
opinion dynamics with an attempt to merge the concepts 
of confidence level and conviction. We find the interest- 
ing result that with a large value of conviction and with 
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any finite bound on the confidence level (i.e., 6 < 1), 
a disordered state exists with tendency to polarisation. 
This is indeed justified, if agents are convinced to a large 
extent in their opinions and interact with like minded 
people only, the sign of the opinion, (representing lik- 
ing/disliking) is likely to be maintained, giving rise to 
a polarised society. The neutral state with all opinions 
equal to zero remains unperturbed with the introduction 
of the bounded confidence. It is found that at least for 
8 > 0.3, the order-disorder transition is first order in na- 
ture. For smaller confidence level, when 6 < 0.3, the 
ordered phase shrinks to a narrow region of the phase 
diagram. 

In conclusion, we obtain a phase diagram with many 
features when the concept of bounded confidence is incor- 



porated in the LCCC model of opinion dynamics. The 
overall result is that when bounded confidence level is 
large, there will be order in the society provided people 
are not too rigid. In the original DNAW model also, it 
had been shown that above a certain confidence level, 
society behaves more homogeneously. We show that this 
tendency remains true but only up to a certain level of 
conviction. This seems to be a realistic scenario and thus 
the combination of concepts from two different models in 
the present model of opinion dynamics is successful in 
reproducing this desired feature of a society. 
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